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Abstract 

We present the high-energy neutrino Monte Carlo event generator ANIS (All Neu- 
trino Interaction Simulation) . The program provides a detailed and flexible neutrino 
event simulation for high-energy neutrino detectors, such as AMANDA, ANTARES 
q^ ■ or ICECUBE. It generates neutrinos of any flavor according to a specified flux 

and propagates them through the Earth. In a final step neutrino interactions are 
simulated within a specified volume. All relevant standard model processes are im- 
plemented. We discuss strengths and limitations of the program. 

o 

""^3 ' Program Summary 

Title of program: ANIS 
1. \ Program is obtainable from: http://www-zeuthen.desy.de/nuastro/anis/anis.html 

Computer on which the program has been thoroughly tested: Intel-Pentium based 
Personal Computers 
Operating system: Linux 
Programming language used: C++ 
^ ■ Memory required to execute: 13 megabyte 

Number of lines in distributed program (version 1.8): 3300 

Libraries used by ANIS: HepMC [1], CLHEP vector package [2] 

Nature of physical problem: Monte Carlo neutrino event generator for high-energy 

neutrino telescopes 

Method of solution: Neutrino events are first sampled according a specified flux, then 

propagated through the Earth and finally are allowed to interact inside a detection 

volume. 

Restrictions of the program: Neutrino energies range from 10 GeV to 10 12 GeV. 

Typical running time: 10 4 events require typically a 1-GHz CPU time of about 300 s 
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1 Introduction 



The first generation of open water/ice Cherenkov neutrino telescopes take 
data (AMANDA [3], BAIKAL [4]) or are under construction (ANTARES [5], 
NESTOR [6]) while the second generation projects such as IceCube [7] are 
already in the planning phase. Such devices aim at detection of cosmic high- 
energy neutrinos. Cherenkov photons from charged particles produced in neu- 
trino interactions near or inside the detector are used to identify neutrinos. 
Because of their large volume and coarse granularity, the energy threshold of 
these open instruments is typically between 10 GeV and 100 GeV and therefore 
high when compared to that of closed underground Cherenkov detectors, such 
as SuperKamiokande and SNO. The flux of neutrinos produced by interactions 
of cosmic rays in the atmosphere (so called atmospheric neutrinos) falls steeply 
with energy (4^- oc E~ 3J ). Hence the event rate due to atmospheric neutrinos 
is largest near the detection threshold energy. This is different for astrophysical 
neutrino fluxes which generally are assumed to have harder spectra. Poten- 
tial sources for astrophysical neutrinos could be astrophysical objects such as 
Active Galactic Nuclei (AGN), Gamma Ray Bursts (GRBs) or heavy relic 
particles, such as produced e.g. in decays of Topological Defects (TD). The 
relevant range of neutrino energies from the more conventional sources such 
as AGN and GRBs range from 10 3 GeV, to 10 9 GeV, energies while the range 
of neutrinos expected from TD models extends up to ~ 10 15 GeV energies. 

Apart from the desired astrophysical neutrinos and atmospheric neutrinos, 
high-energy neutrino telescopes detect muons produced in cosmic ray inter- 
action in the atmosphere (so called atmospheric muons). Due to the low in- 
teraction probability of neutrinos, the trigger rates due to atmospheric muons 
is generally orders of magnitude higher than that of atmospheric neutrinos. 
Hence a large number of search strategies are used, for example using the 
Earth as a shield and selecting only events with horizontal and upward going 
directions. 

In the analysis of the data of neutrino telescopes, Monte Carlo (MC) sim- 
ulations play a decisive role in designing and testing event reconstructions 
and filters, determining filter efficiencies and interpreting the results. Various 
programs exist for event simulations. For example, atmospheric muons can 
conveniently be simulated with the air-shower program CORSIKA [8]. The 
generated muons are further propagated through the overburden consisting of 
water or ice using muon propagation programs such as MMC [9], MUM [10] 
or PROPMU [11]. In a final step the detector response needs to be simulated 
using detector specific programs. Programs for simulating neutrino events are 
known to exist within various collaborations, however, to our knowledge they 
are not publicly available. In particular most programs simulate only some 
specific classes of events (for example v^ reactions) in a limited range of ener- 



gies. 

Here, we present the program ANIS for simulating neutrino events of all fla- 
vors in a wide range of energies. ANIS is a Monte Carlo event generator, 
which prepares neutrinos, propagates them through the Earth and in a last 
step simulates a neutrino interaction within a specified volume around the 
detector. The aim of the program is to provide a tool for precise simulation 
of neutrino events of all flavors in the energy range addressed by high-energy 
neutrino telescopes. The program has been designed to fulfill the requirements 
for analysis of AMANDA data, and could hence be of use for other neutrino 
telescopes. 

An application of ANIS was demonstrated in [12], where the expected rates 
for v e events as calculated with ANIS were presented. The source code as 
well as further documentation can be obtained from the project web site: 
http://www-zeuthen.desy.de/nuastro/anis/anis.html. 

The paper is organized as follows. First a general description of the program 
is given in Section 2. The implemented neutrino interaction processes are 
shortly discussed in Section 3. The possibility to add a new (non-standard 
model) interaction process to ANIS is demonstrated there by an example. 
In Sections 4 and 5 the role of tau leptons and how neutrino propagation 
is simulated in ANIS are discussed. A description of the treatment of the 
final neutrino interactions inside the detection volume follows in Section 6. 
Additional information relevant for potential users is given in the appendix. 



2 Description of the program 



ANIS is written in C++ and uses the vector package of the CLHEP library [2] 
as well as the HepMC Monte Carlo event record [1]. The event record holds 
all interaction vertices with their respective incoming and outgoing particles. 
The currently implemented interaction channels include charged current (CC) 
and neutral current (NC) interactions as well as resonant W~ production: 
u e e~ — > W~ — > anything. 

Primary neutrinos are randomly generated on the surface of the Earth. Cur- 
rently only power-law energy spectra <ft u {E) oc E~ a are implemented, however 
the code is flexible enough to allow its extension to arbitrary neutrino spectra. 

ANIS propagates neutrinos in small steps towards the detector. In interac- 
tions with matter they are either absorbed (CC-case) or regenerated at lower 
energies (NC-case). In the special case of CC v T interaction, a short-living 
r-lepton is produced. It propagates in matter, thereby loosing part of its en- 



ergy, and finally decays giving rise to a secondary v T and, in ~ 17 % of the 
cases, to secondary v^ or v e . The density profile of the Earth used for neutrino 
propagation is chosen according to the Preliminary Earth Model [13]. 

Once the detection volume is reached, a final vertex is generated along the 
neutrino trajectory within the detection volume. In the case of a CC v^N- 
interaction, ANIS correctly simulates the muon scattering angle. 

Along with the full event, four weights are stored: a normalization constant, a 
weight proportional to the total interaction probability of the neutrino as well 
as two weights corresponding to the atmospheric flux of electron and muon 
neutrinos [14]. 

Event rates for atmospheric and various extraterrestrial neutrino spectra are 
obtained by applying the appropriate weights to the events. This last step can 
conveniently be done by a user defined energy dependent weight function ap- 
plied within analysis programs as e.g. PAW [15] and ROOT [16]. Applications 
of the weights are further discussed in appendix D. 

There are currently two output formats implemented: i) the AMANDA specific 
event format J2000 [17] and ii) the HepMC ascii format [1]. 



3 Neutrino cross sections 



At E v < 1 x 10 6 GeV, deep inelastic z/iV-cross-sections may be successfully 
described in the framework of pQCD. Parameterization of the vN structure 
functions, F[ N (x,Q 2 ), may be chosen e.g. according to CTEQ5 [18]. Here, 
—Q 2 = q 2 is the invariant momentum transfer, with the 4-momentum q = 
k u — pi given by the difference in momenta between the incoming neutrino and 
outgoing lepton. The "Bjorken" variable x is given by x = Q 2 /(2Mn(E u — E{)) 
where E v and E\ are the energies of the incoming neutrino and outgoing lepton. 

At higher energies cross-sections are dominated by scattering off sea-quarks 
with small x. The unknown behavior of the structure functions at x < 10~ 6 
makes the calculations model dependent. The uncertainty in extrapolations of 
F" N (x,Q 2 ) to small x and large Q 2 influences the expected cross-sections at 
high energies. Two options are presently included in ANIS: i) the smooth both 
with respect to x and to Q 2 power-law extrapolation of the pQCD CTEQ5 
parameterization to small x and large Q 2 , and ii) hard pomeron [19] enhanced 
extrapolation [20]. The cross-sections of the first case, denoted in the Fig. la 
as pQCD, practically coincide with [21], while the second model (HP, dash- 
dotted curves) predicts approximately 2 times higher cross-sections at E v = 
1 x 10 12 GeV. 
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Fig. 1. a) uN- and resonance v e e -cross-sections and b) normalized ^-distributions. 
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A kinematically important variable is y — 1 
energy transfer to the outgoing lepton and to the final hadronic state. The 
^/-distributions, normalized to the corresponding cross-sections, are plotted in 
Fig. lb; the solid and dashed curves stand for CC and NC z/iV-interactions, 
respectively. The scattering angle, 9, between incoming neutrino and outgoing 
lepton is given by 



cos 9=1 



xy m^f 



i-y E v 



The used cross-section data for CC- and NC-reactions are taken from pre- 
calculated tables. The total cross-section is obtained through interpolation 
between energy bins. A final state can be represented by the pair of x and y. 
Large sets of x and y have been generated for different energies and are stored 
in data tables. During generation of neutrino events such pairs are randomly 
drawn from the tables. The use of pre-calculated tables makes the program 
fast and independent of other packages. 

Since m e <C m^, at high energies practically all z/e _ -cross-sections are neg- 
ligibly small, with the exception of resonant W~ production v e e~ — > W~ — > 
anything at E Pe pa 6.3 • 10 6 GeV (the so called Glashow resonance, [22]) (see 
Fig. la). This process is also included in ANIS. 

ANIS was designed such that new cross-sections and processes can be added 
at a later stage. Assuming for example, one would like to allow microscopic 
black hole production [23] one has to derive a new (C++) class for this process 
and provide the relevant member functions: 



class SigmaBH : public Sigma { 
public: 



double GetSigma( const GenParticle* neutrino, double material []) ; 
void FillVertex(GenVertex* vertex, double material []) ; 

} 

GetSigma returns the cross section, and the function FillVertex inserts the 
final state. By adding this new process to the list of know processes, it will be 
included during simulation. Thereby all particles produced in the final states 
are written to the event record and in case of neutrinos are further propa- 
gated. In this way, any potential implications of non-standard model physics 
can easily be tested. 

Already implemented SM processes can be modified even without modification 
of the code, simply by replacement of the cross-section data tables obtained 
using other programs. 



4 Treatment of tau leptons 



Tau leptons can be produced in CC z/ r iV-reactions and in resonance v e e~ — ► 
W~ — ► v t t~ scattering. The average tau decay length is l r = 49 mx (E T /PeV). 
For energies larger than E T ~ 10 8 GeV, the energy loss during propagation has 
to be taken into account, leading to a slower growth of the decay length. In 
the current version of the program, the stochastic effects in the energy loss 
rate are neglected. The energy loss rate is assumed to be continuous: 

dE 

-^ = a + (3(E)E T , (1) 

with (3 = (1.508 + 6.3 • (E/10 9 GeV) a2 ) • W~ 7 cm 2 /g [24]. 

In contrast to muon propagation, neglecting stochastic energy losses does not 
alter the average r range significantly [25]. At energies of E T ~ 10 8 GeV, the 
continuous case is essentially indistinguishable from the stochastic case and at 
energies around E T ~ 10 11 GeV the difference between the two cases is only 
a few percent. The reason, next to the energetic tau decay, is that in case of 
taus (hard) bremsstrahlung is suppressed while (softer) nuclear interactions 
dominate. 

In ANIS, the tau is propagated in small steps of energy. Once the age in the tau 
rest-frame exceeds the previously determined lifetime the tau decay routines 
are called. 

The tau decay is simulated with the help of TAUOLA [26]. Instead of link- 
ing directly to TAUOLA, we have created tables of final states, consisting 



of the fractional energy: hadronic (tt^jK, ..), electro-magnetic (e, 7r° and 7), 
muons, and the various neutrino flavors. Thereby we have assumed left-handed 
polarization of the tau leptons, which is a good approximation for neutrino 
interactions, if E T /m T ^> 1. The tau decay products are added to the event 
record and in case of neutrinos are further propagated. 



5 Neutrino Propagation 



Neutrinos are generated at the surface of the Earth and then propagated to 
the detector. They are assumed to travel along straight trajectories, neglecting 
the small angle neutrino scattering. This is a good approximation, since at en- 
ergies below 10 5 GeV when the average deflection angle (for example through 
NC reactions) becomes larger 0.1°, the neutrino interaction probability in the 
Earth is small. A deflection angle of 0.1° is below the current instrumental 
resolution of existing or planned neutrino telescopes. The density profile of 
the Earth is chosen according to the Preliminary Earth Model [13]. Neutrinos 
are propagated in small steps through the Earth and in case of an interac- 
tion one of the active interaction channels is chosen according the inclusive 
cross-section. Its final state member function is invoked, which adds the final 
state of the interaction to the event record. Naturally, all propagation effects 
such as neutrino absorption, v T — > t — > v T regeneration, regeneration in NC 
interactions as well as possible neutrino production in the resonant scattering 
of electrons (e.g., V e — > W~ — > v e ,n,r) are included. 



6 Detection volume and final vertex 



The size of the detection volume needs to be specified in the steering file. 
The shape of the volume is a cylinder with its z-axis parallel to the neutrino 
directions. Since open water/ice Cherenkov detectors have a sensitivity to 
events with vertices outside their fiducial volumes, the size of the cylinder 
should be chosen appropriately. In particular due to the potentially long range 
of muon and tau leptons, the final cylinder should extend to the direction of 
the neutrino origin. As an alternative, one can choose a cylinder with a variable 
length which adjusts to the energy-dependent maximally possible muon and 
tau range. The latter allows a more efficient simulation. 

Neutrinos reaching the detection volume are forced to interact within this 
volume. The interaction probability is stored in the event record. Two options 
are implemented: i) all events are written to disk - in order to obtain a physical 
spectrum, one needs to weight all events with their interaction probability, and 
ii) events are sampled according their interaction probability - the events are 
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Fig. 2. Average angle 6 between the incoming neutrino and the outgoing muon in a 
CC reaction as a function of the neutrino energy. 

sampled using an acceptance-rejection method, where the selection criteria are 
determined at run time during initialization. 

The final state depends on the flavor of the interacting neutrino, the type of the 
reaction as well as on the variables x and y of the interaction. Since Cherenkov 
neutrino telescopes cannot distinguish the different types of hadrons, the 
hadronization of the final state is omitted and all hadronic energy is writ- 
ten out as pions. Since neutrino telescopes have a muon angular resolution 
of 0(1°), the scattering angle between incoming neutrino and outgoing muon 
needs to be taken into account. Fig. 2 shows the average scattering angle for 
neutrinos and anti- neutrinos as obtained with ANIS. 

All muons of the final state have to be further propagated using dedicated 
muon propagations programs (such as MMC or MUM), before the events are 
passed to the detector simulation. 



7 Conclusion 



A new program for the simulation of neutrino induced events for high-energy 
neutrino telescopes is presented. The program fulfills current requirements for 
data analysis of neutrino telescopes, and is designed flexible such that new 
processes could be added in the future. The current version of ANIS includes 
all relevant SM neutrino interaction processes. The energies for which cross- 
section data is provided ranges from 10 GeV to 10 12 GeV. 
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A Organization of the Program 



The ANIS source code consists of several files in two directories: anis/ and 
sigma/. The first directory contains the part of the code related to the sim- 
ulation of events, while the latter contains the part of the code related to 
interaction and decay of particles. 

anis/ 

anis.cc - the main file of the program which calls the main processes. 

Steer. cc - the steering file format. 

AnisEvent.cc - the ANIS event record definition. 

FluxDriver.cc - contains classes for simulating neutrino fluxes. 

NeutrTraject.ee - contains classes for calculating neutrino trajectories. 

Propagation.ee - contains classes related to the propagation of neutrinos 

through the Earth to the detector. 

FinalVolume.ee - contains classes related to the final detection volume. 

FinalVertex.ee - classes for sampling the vertex within the final volume. 

Earth. cc - contains information on the Earth density profile. 

AnisIO.cc - definition of various output streams, (see also AnisEvent.cc). 

sigma/ 

Sigma.cc - the virtual base class for all interaction processes. 

SigmaSM.cc - all Standard Model processes. 

SigmaTotal.cc - the container class holding all activated processes. 

ProcessManager.ee - controls the active processes for the neutrino flavors and 

target materials. 

Decaycc - contains all decay classes. 

Table.cc - contains various utility classes for reading the data tables. 



Next we describe the most relevant classes and member functions: 

class FluxDriver 

void Inject (AnisEvent *evt) 

- This method inserts the initial neutrino into the event record. The position 
of the neutrino origin is chosen on the Earth surface, such that the straight 
trajectory of the neutrino intersects the final detection volume. 

class Propagate 

void ThroughEarth (AnisEvent *evt) 

- this method propagates all neutrinos in the event record up to the backside 
of the the final detection volume. 

class SigmaTotal 

double GetSigma (HepMC: :GenParticle *nuin, double material []) 
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void FillVertex (HepMC: :GenVertex *vtx, double material []) 

- The GetSigma member function returns the total cross-section for a given 
neutrino and material (in the current version of the code isospin invariance is 
assumed: all interaction material contains equal numbers of neutrons, protons 
and electrons). The FillVertex member function samples an interaction type 
from the list of possible processes and fills the vertex record with corresponding 
outgoing particles. 

class ParticleDecay 

void AddVertex (HepMC: :GenParticle *parent, AnisEvent *evt) 

- this method inserts a decay vertex for the incoming parent particle into the 
event record. So far the only unstable particle implemented is the tau-lepton. 

class FinalVertex 

void ThrowVertex (AnisEvent *evt) 

- this method inserts a final neutrino interaction (vertex) located within the 
final detection volume into the event record. The interacting neutrino is chosen 
(with a probability proportional to the its interaction probability) from all 
neutrinos reaching the final detection volume. 

Further there is the control class Steer, which holds the run parameter of 
ANIS as well as the ProcessManager which controls the list of available in- 
teraction processes. 

Further documentation of the code can be found on the project web site: 
http://www-zeuthen.desy.de/nuastro/anis/anis.html. 



B The steering file 



The program ANIS reads a steering file called "steer.dat". In case a file name 
is given in the command line, this file is taken as the steering file. The steering 
file contains all run specific information. Here is an example: 



# The flux information 
f p 1 le2 le4 -1. 1. 12 

# The geometry information (volume is cylinder) 
g c 300. 300. 300. 1.7e3 1 

# The run information 
r 100 1 1 2123421 

# The cross-section processes and data. 

s CC cteq5/cteq5_cc_nu.data cteq5/f inal_cteq5_cc_nu.data 101010 110 

s CC cteq5/cteq5_cc_nubar .data cteq5/f inal_cteq5_cc_nubar .data 10101 110 
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s NC cteq5/cteq5_nc_nu.data cteq5/f inal_cteq5_nc_nu.data 101010 110 

s NC cteq5/cteq5_nc_nubar .data cteq5/f inal_cteq5_nc_nubar .data 10101 110 

s GR dummy dummy 10000 1 

# the data directory 

d . . /data 



Lines starting with a letter: s,f,g,r,d specify a program command line. Lines 
starting with # are for comments only. 



f (flux): The f must be followed by a parameter specifying the spectrum type 
(p for power law, this is so far the only implemented one), the spectral 
index, the minimal neutrino energy, the maximal neutrino energy (in GeV), 
the minimal cosine of the zenith angle of the neutrino direction, and the 
maximal allowed cosine of the zenith angle. The last number indicates the 
neutrino type with the following possibilities: 12(= v e + v e ), 14(= z/ M + P M ), 
16(= v r + v T ). 

g (geometry): The g is followed by the type of the final volume (c stands for 
cylinder, which is the only one implemented so far), the radius (in meter) 
and the +/- height of the volume. The height is given by two numbers, the 
first is the positive height (target region), and the second is the height in 
the detector region. The next element is the depth of the detector from the 
surface (in meter). The final number is optional: if set to 1, a geometry op- 
timization is invoked: if the lepton range is larger than the specified volume, 
the volume is adjusted to the maximal lepton range. So far this option is 
implemented only for muon-neutrino. If the number is missing or set to 0, 
then no re-adjustment is made. 

r (run parameter): The r is followed by the number of events to be simulated 
and the output format (l=HepMC,2=f2000). The third value indicates if 
weighted (1) or unweighted (0) events should be simulated. The last number 
is an optional seed for the generator. If no seed is present, the system time 
is used to seed the random number generator. 

s (sigma): The information following the s is the type of the process (CC: 
charged current, NC: neutral current, GR: Glashow resonance) the file con- 
taining the cross-section data, the file containing the final state data, the 
active neutrino bit mask. (y e v e v ll p p y T v T ^ 101010 means all neutrinos, 10101 
means all anti-neutrinos and 10000 means only electron anti neutrinos), 
and the active target material bit mask (110 means protons and neutrons, 
1 means electrons). 

d (data directory) this is the path to the data directory, containing the cross 
section, tau decay and neutrino flux data. 
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Fig. B.l. The geometry: the coordinate system is anchored at the detector center, 
with the z-axis pointing away from the Earth center. The cylinder represents the 
final volume, in which potentially detectable neutrino interactions are simulated. Its 
size is adjustable through the steering file. 



C Output 



Here only the HepMC event output format is discussed. An example for a 
neutrino event as generated with ANIS is shown below. It is the first of 1000 
events generated with the steering file described above. 



GenEvent : #1 ID=1 SignalProcessGenVertex Barcode: 

Current Memory Usage: 1 events, 2 vertices, 3 particles. 

Entries this event: 2 vertices, 3 particles. 

RndmState(0)= 

Wgts(4)=1.20431e-06 5.16006e+17 4.0188e-ll 1.21211e-09 

EventScale -1 GeV alphaQCD=-l alphaQED=-l 

GenParticle Legend 

Barcode PDG ID ( Px, Py, Pz, E ) Stat DecayVtx 

Wertex: -1 ID: 1 (X,cT)=+5.37e+06,+l .26e+06,-2.98e+06,+0.00e+00 

Wgts(l)=0 

0: 1 10001 12 -4.32e+03,-1.01e+03,+2.40e+03,+5.05e+03 1 -2 

Vertex: -2 ID: (X,cT)=-3.50e+02,-7.01e+01,-8.75e+01,+0.00e+00 

I: 1 10001 12 -4.32e+03,-1.01e+03,+2.40e+03,+5.05e+03 1 -2 

0: 2 10002 11 -1.28e+03,-3.00e+02,+7.09e+02,+1.49e+03 1 (nil) 

10003 211 -3.05e+03,-7.15e+02,+1.69e+03,+3.56e+03 1 (nil) 



Events consist of a list of vertices with attached incoming and outgoing par- 
ticles as well as a list of weight. For a description of how to access events, 
vertex, particles and weights we refer to the HepMC documentation [1]. 

• All vertex positions are relative to the detector center (in units of meters). 
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The orientation is such that z axis is pointing upwards. 

The particle identification numbers follow the Particle Data Group scheme. 

Energy and momentum are given in units of GeV. 

Events have four weights attached: the interaction probability in the final 

volume (-Pint), a normalization constant (N) and two weights (i?atm an d 

i?atm) corresponding to the atmospheric neutrino flux for v e and v^ [14]. 



D Using ANIS weights 



To obtain physical distributions of variables (such as neutrino energy or fla- 
vor at the detector) from ANIS events, one can conveniently use the weights 
described in appendix C. 

If events are generated according a spectral index 7 the application of an event 
weight N • Pi n t will result in a simulated flux 



$(£) = (E/GeVHfGeV-V^r-W 



N eve nt 

The normalization is chosen such that N ■ J2 Pint,? results in the total number 

of events per year in the specified volume. Other spectra can be obtained by 
applying additional weight factors. (For example, to obtain a flux with spectral 
index 7', the additional multiplicative weight should be E 1-1 ). 

In case one is interested in an atmospheric neutrino spectrum, each event 
should be weighted according a weight N • P int • Pt atm , where the flux weight 
Patm (which is provided by ANIS) should be chosen according to the simu- 
lated neutrino flavor. Again, the normalization of the resulting distribution 
corresponds to events per year in the specified volume. 

Accounting neutrino oscillations of atmospheric neutrinos can be done using 
weighting functions. Assuming one is interested in muon neutrinos, and P v< ,-> v 
is the energy and direction dependent probability for electron neutrinos to ap- 
pear at the detector site as muon neutrinos, while P Vtl -+ VjX is the probability of 
the initial muon neutrino to be detected as muon neutrino. The corresponding 
weight for each event is then: N ■ P int • (P v v ■ R% m + P v ^ v ■ R% m ). 
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